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By placing an atom into a cavity, the electromagnetic mode structure of the cavity is modified. 
In Cavity QED, one manifestation of this phenomenon is the appearance of a gauge-dependent 
diamagnetic term, known as the contribution. Although in atomic Cavity QED, the resulting 
modification in the eigenmodes is negligible, in recent superconducting circuit realizations, such 
corrections can be observable and may have qualitative implications. We revisit the canonical 
quantization procedure of a circuit QED system consisting of a single superconducting transmon 
qubit coupled to a multimode superconducting microwave resonator. A complete derivation of the 
quantum Hamiltonian of an open circuit QED system consisting of a transmon qubit coupled to 
a leaky transmission line cavity is presented. We introduce a complete set of modes that properly 
conserves the current in the entire structure and present a sum rule for the dipole transition matrix 
elements of a multi-level transmon qubit coupled to a multi-mode cavity. Finally, an effective 
multi-mode Rabi model is derived with coefficients that are given in terms of circuit parameters. 


I. INTRODUCTION 

In single-mode realization of Cavity QED (CQED), a 
single atom coupled to a small high-Q electromagnetic 
resonator can be well-described by a model wherein the 
matter is described by a single atomic transition, and its 
coupling to one of the modes of the resonator can satu¬ 
rate this transition before other modes are populated [1] . 
A plethora of fundamental physical phenomena and their 
recent applications in quantum information science has 
been explored and vigorously pursued with a supercon¬ 
ducting circuit-based realization of this setup [SHI!. In 
such systems, the existence of the atom leads to a modifi¬ 
cation in the cavity modal structure due to Rayleigh-like 
scattering. Such corrections are unobservably small in 
atomic CQED unless a special cavity structure is chosen. 
In recent realizations of circuit quantum electrodynamics 
(cQED) however, such corrections may have observable 
consequences which we discuss in this article. 

A better known manifestation of the aforementioned 
scattering corrections is the so-called term in single¬ 
mode CQED. There has been a lively debate in recent 
years [HHS] about the impact of this term on synthetic 
realizations of the single-mode superradiant phase tran¬ 
sition when instead of one, N identical non¬ 

interacting quantum dipoles are coupled with an iden¬ 
tical strength to a single cavity mode. This particular 
instability of the electromagnetic vacuum has originally 
been discussed [UUTH] within the context of the single¬ 
mode version of the Dicke model m where the term 
was not included. Subsequent work shortly thereafter 
[m? ,iini pointed out that the term rules out such 
a transition. Recent theoretical work on superconduct¬ 
ing realizations of the Dicke Model [8] has challenged 
the validity of such ” no-go” theorems m- Leaving this 
contentious matter aside 0110], we note here that the 
term is a gauge-dependent object, and specifically 
appears in the Coulomb gauge description of the single¬ 
mode atomic CQED. However, the scattering corrections 


due to the existence of an atom in a cavity are physical 
and measurable, and hence not dependent on the choice 
of gauge. In fact, recent realizations m of the mul¬ 
timode strong coupling regime in a very long coplanar 
waveguide cavity, as well as cQED systems in the ultra¬ 
strong coupling regime [221 US] provide settings where 
such corrections may be observable, as we discuss below. 

By fabricating a charge qubit close to a transmission 
line resonator, the resonator’s capacitance per unit length 
is locally altered. This impurity scattering term is typi¬ 
cally neglected [2HlZj in the derivation of the quantized 
Hamiltonian for the multi-mode regime of cQED |2H] , we 
therefore revisit the quantization procedure in section H. 
We discuss how the qubit changes the propagation prop¬ 
erties of the resonator and how as a result this modifies 
its eigenmodes and eigenfrequencies. We show in partic¬ 
ular that this new basis is the one which properly fulfills 
current conservation law at the point of connection to 
the qubit. Finally, in section HI we show how our re¬ 
sults can be generalized to a leaky cavity, one that is 
connected capacitively to external waveguides. In sec¬ 
tion IV, we briefly discuss the comparison to the case 
of atomic CQED. We point out that including the A^ 
term in the Hamiltonian will lead to the same type of 
modification in the modes of a cavity. 


II. MODEL 

As illustrated in Fig. we consider a common cQED 
design |29j consisting of a single transmon qubit that is 
fabricated at point xg (in the dipole approximation, as 
discussed below) inside a superconducting transmission 
resonator of length L. In this section, we assume that 
Cl.b — 0, which corresponds to closed (perfectly reflect¬ 
ing) boundary conditions at a:; = 0, L. In section HI we 
discuss the open case, where a finite transmission line of 
length L is coupled through nonzero capacitors Cji^L to 
the rest of the circuit. 
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FIG. 1. A superconducting transmon qubit coupled capaci- 
tively to a superconducting coplanar transmission line. (a). 
Device view (b). Equivalent circuit 


A. Classical Hamiltonian and CC Basis 


As discussed in detail in Appendix B, the Hamiltonian 
for this circuit can be written as 


Si("'If) 

^^-^ 
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The notation used here follows the canonical approach 
to quantization of superconducting electrical circuits 
[27ll30j . briefly reviewed for completeness in Appendix A. 
In the above Hamiltonian, the canonical variables <i>j and 
Qj represent the flux and charge of the transmon qubit 
respectively. In a similar manner, the canonical fields 
$(a:, t) and p{x, t) represent the flux field and charge den¬ 
sity field of the transmission line. Furthermore, (j)Q = ^ 

is the flux quantum and 7 = ^ . 

There is a crucial difference between the Hamiltonian 
we have found here, with respect to earlier treatments 
[m US]. We do include the modification in the res¬ 
onator’s capacitance per length at the qubit connection 
point xq, c(x, xO) = c + CeS(x — xq ) where is the se- 
ries capacitance of Cj and Cg given as ■ The Dirac 


delta function is the result of treating the qubit as a point 
object with respect to the resonator, whereas a more re¬ 
alistic model would replace that with a smooth function 
as we have discussed in Appendix F in great detail. As 
we see shortly, the delta function appearing in the de¬ 
nominator will not cause any issues in the quantization 
procedure, since the charge density p{x, t) also contains 
the appropriate information regarding this point object 
so that turns out to be a continuous function in 

c(x,xo) 

X. Once we understand how this correction influences the 
photonic mode structure of the resonator, we will move 
on to use that information in constructing the quantiza¬ 
tion procedure. 

The Hamiltonian equations of motion derived from the 
Hamiltonian of the modified resonator including 

the impurity scattering term is given by (see Appendix 
B): 

dt c{x,xo) 

dp{x,t) ^ 1 d'^<^{x,t) 

dt I dx^ ^ ’ 

The solution to these linear equations can be writ¬ 
ten in terms of the Fourier transform ^{x,t) = 

^ d,te~^‘^*^{x,ijj) where $(a;,w) is the solution of 

the ID Helmholtz equation 

+ lc{x,xo)uj'^^{x,u}) = 0 (4) 

We look for solutions that carry zero current across 
the boundaries, implemented by Neumann-type bound¬ 
ary conditions dx^{x)\x=o,L = 0. A solution then exists 
only at discrete and real values w = a;„. The Dirac delta 
function hidden in c{x, Xq) can be translated into discon¬ 
tinuity in dx^ix) which is proportional to the current 

I{x) = — j that enters and exits the point of con¬ 
nection to the transmon 


1 9$(a:,a;) 1 9l>(x,w) 

I dx I dx 

x+ 


CsU}‘^^{xo,uj) (5) 


where the R.H.S is the current that enters Cj through 
Cg , therefore the series capacitance Cg- This condi¬ 
tion amounts to the conservation of current at the point 
of connection to the qubit and thus it is appropriate 
to call the set of eigenmodes satisfying this condition 
the current-conserving (CC) basis. The solution of the 
above-stated Neumann problem gives the CC eigenfre- 
quencies a;„ through the transcendental equation 


sin(fc„L) -I- XsknL cos(fc„xo) cos(A:„(L - xq)) = 0 ( 6 ) 


In the equation above, knL = = '/icoJnL is the nor¬ 
malized eigenfrequency and ^ is a unitless mea¬ 

sure of the transmon-induced modification in eigenfre- 
quencies/eigenstates compared to the conventional cosine 
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FIG. 2. a) The first 10 modified resonances for different values 
of Xs and xq = O.OIL. Higher modes with larger Xs experience 
a larger shift in frequency. b,c) Normalized level-spacing for 
Xs = 0.001, and 0.01 respectively. The blue dashed line shows 
the constant level spacing for Xs = 0. 



FIG. 3. Normalized energy density of the first 4 modes for 
xo = 0.011/ and x® = 0.1. The black curve shows cosine 
modes while the red curve represents GG modes. The blue 
star shows where the qubit is connected. 


basis. The CC eigenfunctions are given by 


cos (kn{L — cco)) cos (knx) 
cos (knXo) cos (kn(L — x)) 


0 < X < Xq 
Xq < X < L 


(7) 


The proportionality constant has to be set by the or¬ 
thogonality relations that can be found directly from the 
modified wave equation as 



^^ c{x, Jo) ^ 



d^mjx) d^njx) 
dx dx 


kmkriLSn 


( 8 ) 

(9) 


Based on these results, eigenfrequencies are not only 
sensitive to Xs, but also to the point of connection xq. 
In order to understand this modification better, first we 
have plotted the normalized eigenfrequencies in Fig. 
for different values of Xs and the case where qubit is con¬ 
nected very closely to one of the ends. This is a standard 
location for fabricating a qubit [21] to attain a strong 
coupling strength between the resonator modes and the 
qubit, since the electromagnetic energy concentration is 
generally highest near the ends. In this figure, the blue 
circles representing the eigenfrequencies for Xs = 0 are 
located at rnr. For Xs ^ 0, all lower CC eigenfrequencies 
are red-shifted with respect to the Xs = 0 solutions and 
by going to higher mode number and higher Xs, the de¬ 
viation becomes more visible. In a larger scale however, 
the behavior of Xs 0 eigenvalues are non-monotonic 
and most notably, display a dispersion in frequency. For 


better visibility of this periodic behavior, in Figs 2bpc 


we have compared the level spacing of CC modes for dif¬ 
ferent values of Xs to the constant level spacing of un¬ 
modified cosine modes. This behavior is determined by 
the position of the qubit connection point Xq and is easy 
to understand. Since xg = O.OIL sits at the local minima 
of modes 50, 150, 250 and so on, we expect a periodic be¬ 
havior in the values of CC eigenfrequencies where within 
some portion of that period set by Xs, CC solutions are 
less than the Xs = 0 solutions and vice versa in the re¬ 
maining portion. 


In Fig. I^we show the spatial dependence of the first 
four modes. The amplitude of the lower CC eigenmodes 
at the qubit location are consistently less than that for 
the unmodified cosine eigenmodes, suggesting that the 
actual coupling strengths of the qubit to these modes are 
below the ones predicted by the Xs = 0 modes. 


We note that similar modification in the modal struc¬ 
ture of resonators have also been noted for transmission¬ 
line resonators containing inline transmons |31j as well 
as 3D cQED architectures |32]. In these studies, the 
modifications result from solving the spectral problem 
of the quadratic Hamiltonian that in addition to the 
resonator part also includes the linear part of the JJ’s 
cos(27r$j/(/)o)-type non-linearity, unlike the situation de¬ 
scribed here. 


In order to see the dependence of these CC modes on 
Xq, we have also investigated two different cases of con¬ 
necting the qubit to Xg = j and a;o = ^ Appendix C 
in figures and respectively. 
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B. Canonical Quantization 


can be expressed as 


As discussed in Appendix C, the conjugate quantum 
fields of the resonator can be expanded in terms of CC 
basis as 

Hx, t) = Yl ( 217 ^) 

p{x,t) = -iY^ (^) c{x,XQ)^n{x) 

( 11 ) 


7^ — ^ A ^ ) hujyiU^Q,n 


n n 



m<n^l 


where the gmni stands for the coupling strength between 
mode I of the resonator and the transition dipole Pmn 
and is obtained as 


By substituting these expressions into 77™°“^ and em¬ 
ploying orthogonality conditions and one finds its 
diagonal representation 

= + ( 12 ) 

T3. X 


In a similar manner, the qubit flux operator c&j and 
charge operator Qj can be represented in terms of eigen- 
modes of transmon Hamiltonian as 

= E ("^1 \n) Pmn{t) (13) 

m,n 

Qjit) =^{MQjW\n) Pmnit) (14) 

m,n 

where Pmnp) represent a set of projection operators act¬ 
ing between states \m) and |n). Working in the flux basis, 
the eigenmodes are found through solving a Schrddinger 
equation as 


' 2Cj 


— Ej cos ( 27r^ 


4>o /. 


(15) 


^9mnl — q 



PQj, mn )‘^l{xo) 


(19) 


Various TRK sum rules [33] can be developed for a 
transmon qubit, as discussed in detail in Appendix D. 
For instance, the sum of transition matrix elements of 
Qj between the ground state and all the excited states 
obey 


E 2 (£;, 

n>0 


l^o) ( 0 |Qj|n) 


{2epEj (0| cos 



| 0 ) 


< {2epEj 


( 20 ) 


Since all terms on the R.H.S are positive, this imposes 
an upper bound to the strength of Qjpn- A multi-mode 
Rabi Hamiltonian can be recovered by truncating the 
transition matrix elements to only one relevant quasi¬ 
resonant transition term (assumed here to be the 0 —?► 1 
transition): 


H = -fiwoio-"' -t- E f^ndhdn 

" ( 21 ) 

+ E + dl) 

n 


whose solution can be characterized in terms of Mathieu 
functions [29]. Due to the invariance of the Hamiltonian 
under flux parity transformation, the eigenmodes are ei¬ 
ther even or odd functions of and as explained in 
Appendix C, only off-diagonal elements between states 
having different parities are non-zero. Consequently, flux 
and charge matrix elements are purely real and imagi¬ 
nary. Therefore, we can rewrite 

h(t) = E (™l ^J(O) 1 ^) (^mn(t) + Pump)) (16) 

m<.n 

QjP) = E (™l Qj(0) (^PmnP) - P„m(0) (17) 

m<.n 


Finally, applying the unitary transformation a; —>■ idi 
and Pmn iPmn for m < n, the second quantized rep¬ 
resentation of the Hamiltonian in its most general form 


The coupling strength is now reduced to 


/ \ 2 

%n=7f^j {iQjfil)^n{xo) 


and based on Eq. 20 Qyoi has to satisfy 
2e^Ej 2Ej 


IQj,oi| < 


El — Eq pSEjEc — Ec 


( 22 ) 


(23) 


2 

where we have defined the charging energy Ec = 

In order to understand how much gn can deviate in 
practice from its former widely used expression in terms 
of the unmodified {xs = 0) modes, in Fig. we have 
compared the results for various values of Xs- We note 
that in recent experiments on an ultra-long (^ 70 cm) 
transmission line cavity m, Xs was found to be around 
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at the beginning of each cycle (black, blue and purple) 
rather than in the middle (red and blue). 

It remains to be seen whether the non-zero dispersion 
in frequencies and the modifications of coupling strengths 
to higher order modes is observable in practice, because 
in considering such large frequency intervals, the fre¬ 
quency dependent response function of superconductors 
would have to be taken into account [HES]. 

The dependence of on xq has been studied in Ap¬ 
pendix C for two other cases xq = ^ and xq = ^ along 
with their corresponding CC eigenfrequencies and eigen- 
modes. 
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III. GENERALIZATION TO AN OPEN-CAVITY: 

OPEN-BOUNDARY CC BASIS 

We now discuss the quantization in an open geometry, 
where the resonator is coupled to two long microwave 
transmission lines, of length and L^, at each side 
through nonzero capacitors Cl and Cn (Fig. [^. In Ap¬ 
pendix E, we discuss how these nonzero capacitances al¬ 
ter the boundary conditions at each end, and hence the 
mode structure as a result as well. The resulting real 
eigenfrequencies of the resonator can be found from the 
transcendental equation 


FIG. 4. Normalized coupling strength for xq = O.OIL 
a) First 20 modes b) Large scale behavior for 250 modes. 
In both graphs, coupling strength is normalized such that 
only the normalized photonic dependence is kept i.e. Qn = 
(^nL) ^ (:ro) 


10“^. For shorter, more standard transmission line res¬ 
onators we should expect Xs ~ O.I because Xs oc 

As we observe in these figures, CC couplings are 
very sensitive to a change in Xs- For instance, in Fig. 
which is for the common case of connecting the qubit 
to one end, xg = O.OIL, it is observed that even for 
Xs = 0.001 (red stars) the relative shift in the highest 
mode shown (mode 20) is about 3%. This relative change 
increases to 26%, 69% and 80% for Xs = 0.01, Xs = 0.05 
and Xs =0.1 respectively. 

These modifications, even for small Xst are clearly ob¬ 
servable in the multi-mode regime, i.e. when the qubit is 
resonant with a very high order mode. To study the large 
scale behavior of couplings, we have plotted the first 250 
CC couplings gn in Fig. |^for the same parameters. As 
we mentioned earlier, due to the fact that the qubit is 
placed at a symmetry point, we expect that with a pe¬ 
riod of 100 modes, the couplings fall down to zero. The 
first mode that has a local minimum at Xg = O.OIL is 
mode 50 and it occurs again at modes 150, 250 and so 
on. As a general rule, higher CC modes experience a big¬ 
ger shift in their coupling strength. Another important 
observation is the suppression of coupling strength as Xs 
increases such that the highest coupling strengths occur 


+ (1 - XRXhiknLf) sin (fc„L) 

+ iXR + Xl) knL cos {knL) 

-I- Xsk„L cos (knXg) COS {kn{L - Xg)) 

(24) 

- XRXsiknL)'^ COS (knXg) sm (fc„(L - Xo)) 

- XLXsiknL)'^ sin (k^Xg) cos (/c„(L - Xg)) 

+ XRXLXsiknL)^ sin (knXg) sin {kn{L - Xg)) = 0 

where xr,l = arc normalized coupling constants to 
the left and right transmission line. Considering only the 
first two terms in the expression above, by setting Xs = 0, 
would lead to the well-known equation in the literature 
[21126] 


tan (knL) 


{XR + XL)knL 

XRXLiknL)'^ - 1 


(25) 


which only describes eigenfrequencies of an isolated res¬ 
onator and does not contain appropriate current conser¬ 
vation at the qubit location. The third term is the same 
modification we have found in the closed case and has 
a significant influence as Xs increases, while the others 
represent higher order corrections and are almost negli¬ 
gible except for very high order modes. The real-space 
representation of these eigenmodes are found as 


$„(a;) oc 


^n{x) 


0 < X < Xg 
Xg < X < L 


(26) 
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where and $>( 0 ;) are given by 

$<(a;) = [cos {kn{L - xo)) - Xi^fc„Lsin (fc„(L - a;o))] 

X [cos (knx) - XLknL sin (knx)] 

(27) 

$>(a;) = [cos(fc„a:o) - xlA:„L sin (fc„a;o)] 

X [cos (kniL - x)) - XRknL sin (A:„(L - x))] 

(28) 


The open-boundary CC basis can be shown to satisfy the 
modified orthogonality relations 


^ C-op{x^ Xq) 


$m(a:)$„(x) = 


(29) 


where the capacitance per unit length Cop(x,xo), due to 
the leaky boundary is given by 


Cop(x,Xo) = c + Cs5{x - Xo) 


^ Cr5{x - L-) + Cl5{x - Q+) 


(30) 


The remaining orthogonality relations for the current is 
also modified, 




i9x 


dx 


- 2 + *n) L [xfll>™(L-)l>„(L-) + xlI>™(0+)I>„(0+) 

— kmknd-'^mn 

(31) 


The same argument holds for the CC modes of the left 
and right transmission lines, while the exact knowledge 
of these modes requires assigning appropriate boundary 
conditions at their outer boundaries. For instance, if the 
side resonators are assumed to be very long, an outgo¬ 
ing boundary condition is a very good approximation, 
since the time scale by which the escaped signal bounces 
back and reaches the original resonator is much larger 
than the round-trip time of the central resonator. On 
the other hand, if we have a lattice [24l - f26| of identi¬ 
cal resonators each connected to a qubit and capaci- 
tively coupled to each other, then the same basis can 
be used for each of them. Assuming we also have the 
solution for the CC basis of right and left resonators as 
e N, S' = {i?, L}}, the quantum flux fields 
in each side resonator can be expanded in terms of these 
CC modes as 

1 

*»(^-‘)=i:( ( )'(vs«)+ic«))Axx) 

n,S ^ ' 

(32) 


where &r!,,s={fi,L} sere the annihilation and creation oper¬ 
ators for the open CC mode in each side resonator. 
Following the quantization procedure discussed in Ap¬ 
pendix E, we find the Hamiltonian in its 2nd quantized 



FIG. 5. A single-electron atom interacting with the EM field 
inside a closed cavity of length L. 


representation as 

^ n _^ ^ n _^ n,S—{L,R} 

■Ha 

T ^ ^Prnn d” “f ^ 

^ "V ^ 

Uirrt 

^ ^ ^Pmn,S (^m “1“ (j^n,S ^n,S^ 

m,n,S—{L,R} 

^-V-^ 

Rcb 

(33) 


In this expression, (3rnn,s stand for coupling strength 
of nrd^ open CC mode of the resonator to open CC 
mode of the side baths and is found as 


Cs 


/3mn,S = )$„,s(L+) 


(34) 


where Cs here stands for side capacitors Cr^l and 
shouldn’t be confused with the series capacitance Cs 
introduced earlier. Notice that light-matter coupling 
strength gmni has the same form as before, but in terms 
of open CC eigenmodes and eigenfrequencies. 


IV. DISCUSSION 

The corrections to the spectral structure of the res¬ 
onator found in sections II and III are mathemati¬ 
cally equivalent to the scattering corrections that re¬ 
sult from the presence of an atom in atomic CQED sys¬ 
tems. Electromagnetic field quantization has been stud¬ 
ied in great detail for CQED systems including single¬ 
electron atoms [31 [37], multi-electron atoms [33], and 
for atoms embedded in dispersive and absorptive di¬ 
electric media [38ll43j . For completeness, in Appendix 
G we present a full derivation of the minimal coupling 
Hamiltonian(neglecting electron’s spin) for this system 
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starting from a Lagrangian formalism that yields the 
Maxwell’s equations and the Lorentz force law [38] . The 

term - $(a;o, t)) that ap¬ 

pears in the canonical quantization of cQED systems 
is mathematically equivalent to the approximate(zero- 
order dipole approximated) minimal coupling term Te ~ 
2 ^{Pe — eA(Ro,t))^ appearing in CQED Hamiltonian, 
thus their impact on the cavity modal structure is sim¬ 
ilar. It could be argued that the freedom in the choice 
of the point of reference for the generalized fluxes i.e. 
the choice of ground, is analogous to the gauge freedom. 
However, the fact that the cavity modes are modified due 
to the existence of the qubit is a property that is gauge- 
independent. In Appendix G, we show that in a similar 
manner to the discussion here, the existence of the 
term in the Coulomb gauge gives rise to modified spectral 
properties of the resonator. However, in atomic CQED 
these corrections are tiny because of the smallness of typ¬ 
ical atomic transition dipoles and the fine structure con¬ 
stant. In Appendix B.3, we have also studied the reverse 
question and proved it is feasible to retrieve an A^-like 
term if one naively performs the quantization by the co¬ 
sine basis. This completes the similarity between cQED 
and CQED Hamiltonians in the lowest order(zeroth order 
dipole approximation) where the dimension of transmon 
(atom) is completely neglected compared to the cavity’s 
wavelength. 


a certain loop, the algebraic sum of flux variables over 
that loop should be equal to the external flux. Taking 
into account all these considerations, the Lagrangian for 
any cQED system is found as 

£[$„,$„] (Al) 

where T represents the kinetic energy corresponding to 
capacitors as Tc[^\ = ^ stands for the poten¬ 

tial energy corresponding to inductors as Ul{^} = 
or any other nonlinear magnetic device such as Joseph- 

son junction Uj[^j] = —Ej cos where <^o = ^ 

is the flux quantum. 


Appendix B: CLASSICAL HAMILTONIAN AND 
MODIFIED EIGENMODES OF A CLOSED cQED 
SYSTEM 

Here, we follow the procedure discussed in Appendix 
A for the system shown in Fig. |lb| W e first use a dis¬ 
cretized lumped element LC-model[44j for the microwave 
resonator and then take the limit where these infinites- 
simal elements go to zero while leaving the capacitance 
and inductance per length of the resonator invariant. 
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Appendix A: cQED NOTATION 

In order to describe the dynamics of any cQED sys¬ 
tem we follow a common quantization procedure HZIED]. 
The first step is to write the Lagrangian in terms of a 
generalized coordinate £[(/„]. Then, a Legendre trans¬ 
formation to find the Hamiltonian in terms of the co¬ 
ordinate and its conjugate momentum as 

TLlgn,Pn\ = X) gnPn ~ ^ ■ Finally, we need to apply the 

n 

canonical quantization by imposing a nonzero commuta¬ 
tion relation between the conjugate pairs as [qn,Pn\ = 
Here we go after the convention used in cQED by choos¬ 
ing the generalized coordinate as <i>„(t) = 
in which 14 , (t) is the voltage at node n and is measured 
with respect to a ground node. This quantity has the 
units of magnetic flux and it can be shown that its con¬ 
jugate variable has the units of charge and we denote it 
by Qn{t)- There is an additional rule one has to keep in 
mind. In the case of external magnetic flux applied on 


1. Discrete Limit 


a. Classical Lagrangian for the Discretized Circuit 


In terms of the generalized coordinates introduced in 
Appendix A, the Lagrangian for the discretized circuit 
can be written as the difference between kinetic capaci¬ 
tive energy and potential inductive energy and it reads 


-b Ej cos ( 27r^ ) 


00 ) 


E 


1 A T2 1 

-cAx^z. — — 

2 2lAx 


($„+i - 4>„)" 


(HI) 


+ - ®j)^ 


In the expression above, we have labeled the discrete 
nodes such that the qubit is connected to the zeroth node. 


b. Classical Hamiltonian for the Discretized Circuit 

The first step to find the Hamiltonian is to derive 
the conjugate variables associated with the generalized 
coordinate {{$„};$,/}. These conjugate variables will 
have the dimension of charge and we represent them 
as {{Qn}',Qj}- By definition, these conjugate variables 





read 


5L 


Qj=-^ = {Cj + Cg)^j - Y, 


5$ 

5L 


(B2) 


Qn=^ = (cAx + Cg6no)^n " CgSnO^J (B3) 

O^ri 

The next step is to calculate the discrete Hamiltonian 
by a Legendre transformation 


H = YQn^n+QAj-C 


— 2 ^9)^J ~ <^ 0 ® ( 27’' — 


E 


-{cAx + CgSno)^n + 


1 


21 Ax 


($„+i - $„)" 




(B4) 


Now, we need to solve for $./ and $„ in terms of Qj 
and Qn to represent the Hamiltonian only in terms of 
generalized coordinates and their conjugate variables. 

Before proceeding further, let’s define a few quantities 
in order to simplify the calculation 




7 = 


Cg+Cj 


c,= 


CgCj 

Cg+Cj 


Cs^n — cAx -\- Cs^nO 
^g,n = cAx + CgSnO 


(B5) 

(B6) 

(B7) 

(B8) 


where Cs represents the series combination of the cou¬ 
pling capacitor Cg and Transmon’s capacitor Cj. In 
terms of these new quantities we can write 


= 


Qn 'y^noQj 


C: 


+ 


s,n 

1 

Cj 


a 




s,0 


iSnO 

r' 

. rj. 


Qn 




n = COS {27Ti^ 


2 7 


$0 


E 


2^9,nK 


1 


21 Ax 


($„+i - $„) 


(B9) 

(BIO) 

(Bll) 


- Cg^O^J 


By inserting the expressions for <i)„ and into the 
one for the Hamiltonian one finds that 




n 

_7 

C«.o 


J_ , 7 ^(Cg.O - iCg) 

' r^2 

^9 ^s ,0 

1 C'g,« - iCgSnO ^2 

2 

QjQo 


-Qt. 


9 / 4>7 

Qj - Ej cos ( 27r— 
2lAx 


(B12) 


8 

Notice that this expression can be further simplihed 
since Cg^n and Cs,n are related as Cg^n — iCgSno = Cs^n- 
Therefore, the final result for the discretized Hamiltonian 
will be 





(5 j - Ej cos 





1 

21 Ax 


(4^n+l 



+ -^—QjQo 


(B13) 


where the conjugate variables obey the classical Poisson- 
bracket relations 


\^mQm} — ^mn (B14) 

{4>j,gj} = l (B15) 

{Qn,Qm} = {^n,^m }=0 (B16) 

{gj,gj} = {4>j,$j} = o (B17) 


2. Continuum Limit 


Now that we have the expressions for Lagrangian and 
Hamiltonian in the discrete limit, we can obtain the 
analogous continuous ones by simply taking the limit 
Ax —>■ 0, while keeping the capacitance and inductance 
per length constant. In order to do so, let’s find first how 
some of the terms change in this limit. Let’s first inves¬ 
tigate the Kronecker delta. It is quite natural to argue 
that 


I. ^nO 

lim —— 

Ax-j-o Aa; 


(5(a;) 


(B18) 


6{x) here represents the Dirac delta function. One can 
simply verify that it has all the properties of a Dirac delta 
function 


1. S{x) = 0, a: 7 ^ 0 

2 . 6{x) —>■ -boo, a: —>■ 0 

3. lim S(x) dx = 1 


Based on this result, it is possible to find how Cg^n 
transform in the continuous case as 

c(a;) = lim = c + Cs5{x) (B19) 
As-s-O Ax 


We call this quantity modified capacitance per length of 
the resonator, since it has the information regarding the 
position of the qubit and the way it changes the capaci¬ 
tance at the point of connection. By going to continuum 
limit, the charge variable Qn goes to zero, since it rep¬ 
resents the charge of infinitesimal capacitors. However, 
the charge density remains a finite quantity 


p{x,t) 


lim 
Arc—>-0 


Qnjt) 

Ax 


(B20) 
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Finally, by definition 

d>„+i(t) - $„(t) 


lim 

tc—>-0 


Ax 


9$(x, t) 
dx 


(B21) 


a. Classical Lagrangian and Euler-Lagrange E.O.M. in the 
Continuum Limit 


In terms of voltage and current, equa tion |B29 means 
that voltage is spatially continuous while B30| means that 
current is not continuous at the position of the transmon, 
since some of the current has to go into the qubit. The 
two terms on the right hand side of |B30| are proportional 
to the current that enters Cj and the Josephson junction 
respectively. 


Applying the limits introduced in the previous section, 
Lagrangian in the continuum limit reads 


h. Classical Hamiltonian and Heisenberg E.O.M. in the 
Continuum Limit 


C=l{Cj + Cg)^^J-Uji^j) 


.L/2 

+ / dx 

J-L/2 

— dx CgS{x)^j 

J-L/2 


1 - ^ r/ 1 

2(c + Q«(x))(-) 


dt 


(B22) 


where Uj{^,j) = —Ej cos . Euler-Lagrange 

equations of motion are derived from the variational prin¬ 
ciple = 0 as 


iCg + Cg)^g - Cg dxSix)^ + = 0 

(B23) 




$.7 =0 


(B24) 


It is h elpfu l to rewrite these equat ions by first finding 


from B23 and plugging into 


B24 




(B25) 


which is a wave equation with modified capacitance per 
length and the transmon qubit as a source on the right 
hand side. Therefore, the simplified equations of motion 
read 


- , 7 dUj{<^j) 52$(0,t) 

^ Hi - ^ - = 7 


Cg 9 $, 


dE 

dU.j{^.j) 


(B26) 

(B27) 


The Dirac delta function in the wave equation |B27 


can 


be translated into discontinuity in the spatial derivative 
of 4>(x, t). Therefore, equation B27 can be understood as 


52 $ cl2(j) 

^-ic^=0,x,i0 

$(0+,t) = $(0-,t) 


dx 


c=0+ 


dx 


= IC.. 


x—0~ 


92 $ 

dE 


■h 


x—0 


(B28) 

(B29) 

dUji^.j) 

d^j 

(B30) 


Starting from our discrete Hamiltonian, we try to find 
the continuous one again by taking the limit Ax —> 0. 
Let’s investigate each term seperately. 


lim ( -Z- 


Ax—^O \ C, 


1 

Cs,0 


1 


Therefore, the transmon’s Hamiltonian will be 
-Ej cos (2 tt^ 


SCj V J’o. 

The resonator’s Hamiltonian transforms as 

'Ql 1 


(B31) 


(B32) 


lim y 

V on — ‘ ^ 


Ate—>-0 




2 c. 


= lim > Ax 

A T' _ 


Aai—^0 
rLl2 


2lAx 

Ks) 


dx 


Ax 

p^{x,t) 1 


1 

Jl 


$n+l - 


d^{x, t) 


Ax 

2" 


(B33) 


-L/2 2c(x) 21 V dx 

and finally the interaction term can be written as 
7 ^ ^ Qr. 


t^QjQo = i™„7Qj y 

Ax—>-0 C/g Q Ax—>-0 


Cs 


JnO 


= lim 7 Q 7 > 
Ax^O ^ ^ 


Qr^ A 
Ax ""0 


Ax 

n Ax 

=^0, A" iCALsi,) 

J-L/2 


Ax (B34) 


c(x) 


Putting all the terms together, the final expression for 
the Hamiltonian will be 


7 L /2 


dx 


-L/2 


Ha 

2c(x) 


1 / 94>(x, t) 
21 V dx 


(B35) 




7Qj dx^^6{x) 
J-L /2 C(x) 

-V-" 

H,„t 
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where the Poisson bracket relations now change to 

= l (B36) 

p{x',t)} = S{x — x') (B37) 

Notice that in our expression for Hamiltonian we have 
a Dirac delta function hidden in c{x) in the denominator 
of both resonator’s capacitive energy and the interaction 
term. At the first sight, it might seem unconventional to 
have a Dirac delta function in the denominator. However, 
we will show that the charge density p(x, t) is also pro¬ 
portional to c(x) which makes these integrals have finite 
values. 

we know that the time dependence of an operator 
o ({^’n}, {Qn}; ^"7, Qj; t) is determined by 


By combining the above equations and rewriting 
them in Fourier representation in terms of <I>(x,w) = 
dt^{x, we obtain 

-I-k(x,xo)w^l>(x,a;) = 0 (B46) 

Notice that there is a Dirac delta function hidden in 
c(x,Xo). As we mentioned earlier, this can be translated 
into discontinuity in dx^{x) which is proportional to the 

current I{x) = that enters and exits the point 

of connection to the transmon 


1 9I>(x,Cli) 1 d^{x,uj) 

I dx I dx 

x'X 


C'sClI^$(xo,w) 

(B47) 


Using the Poisson-bracket relations introduced above 
one can find the Hamiltonian E.O.M as follows 


94>(x, t) 
Wt 

dp{x,t) 
dt 
94>j 
dt 

dQj 
dt 


p{x,t) ^ i5{x)^ 
■ —rx—I— 
c{x) c{x) 

1 t) 


l dx"^ 


= ^ + r 


Cj J-L/2 

c(x) 

dUj{<i>j) 



00 


p{x,t) 


(B39) 

(B40) 

(B41) 


We are after a complete set of modes $n(x) = <i)(x, a;„) 
where any solutions to the previous wave equation can 
be linearly decomposed on them. In order to find these 
modes, we have to solve 


d^^n{x) 


dx"^ 

with boundary conditions 


-I- Icoj^^nix) =0, X ^ Xo 



(B42) 

d^n{x) 


d^n{x) 

V 00 / 


dx 

X—Q 

dx 


= 0 


(B48) 


(B49) 


c—L 


The results here, can be generalized to a case where 
the transmon is connected to some arbitrary point xg, 
where the modified capacitance per length now changes 
to c(x, Xg) = C + Cs6{x — Xg). 


d^n(x) 

dx 


d^n(x) 

dx 


+ lCsUji^n{xg) 




0 (B50) 
(B51) 


3. Modified Resonator Eigenmodes and 
Eigenfrequencies 


Consider the second term in B35| which is the 


modified resonator Hamiltonian. The goal here is to find 
out how this modification in capacitance per length in¬ 
fluences the closed Hermitian eigenmodes and eigenfre¬ 
quencies of the resonator. Assuming that the transmon 
is connected to some arbitrary point Xg the Hamiltonian 
is given as 


Applying the boundary conditions we find a tran¬ 
scendental equation whose roots will give the Hermitian 
eigenfrequencies of this closed system as 

sin(fc„L) -I- XskjiL cos{knXg) cos(A:„(L — Xg)) = 0 (B52) 

In the expression above, kn represents the wavevector 
defined as k^ = Iclo^ and the quantity Xs = ^ is 
a unitless measure for the discontinuity of current 
introduced by the transmon. 


'TjTnod 

He 


dL 

1 dx 

p^{x,t) ^ 1 

fd^x,t)Y 

h 

2c(x,Xg) 21 

V ax ; J 


(B43) 


Applying the Poisson-braket relations discussed in the 
previous section, the Hamiltonian E.O.M for the conju¬ 
gate fields read 


pix,t) 
dt c(x,xg) 

dp(x,t) ld^<l>(x,t) 

dt 1 dx^ 


(B44) 

(B45) 


Eventually, the eigenfunctions are found as 


cos {kn{L — Xg)) COS (fc„x) 
cos (knXg) COS (/Cn(T — x)) 


0 < X < Xg 
Xg < X < L 

(B53) 


where the proportionality constant is set by the orthog¬ 
onality relation 

(B54) 

C 
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Another important orthogonality condition can be de¬ 
rived in terms between {dx^n} as 


d^Md^r^_k k L5 

LiUy ^ ^ — rvjYirxjn-^'^mn 


dx 


dx 


(B55) 


Finally, it is instructive to show explicitly the origin of 
an A^-like term when instead of the CC-basis the con¬ 
ventional cosine modes are chosen. Replacing p{x, t) from 


B44 


in 


nimod 

He — 


gjygg 

\c{x,Xo) (d^{x,t) 




21 


dx 


(B56) 


erators and p{x,t) to be 

1 

h 


Hx,t) = 

n 

p{x,t) = -iY^ 


2uJnCL 

f hxtn 
\2cL 


(a„ + aj,) l>„{x) (C5) 

(a„ - aD c{x, a;o)l>„(x) (C6) 


Substituting these expressions into Wq°'^ and using the 
orthogonality relations B54 and B55 will result in 


n 


mod 


huj> 


= X! = X! Indian + const. 

n 71 

(C7) 


Substituting c(x, Xg) = c + CsS{x — xq) leads to 


2^ di~ J dx J 

'' -,,-' 

He 

, ( 5$(a;o,t)V 

(B57) 



which is a sum over energy of each independent mode as 
we expected. Having found the the resonator’s Hamilto¬ 
nian in second quantized form, we have to calculate now 
the spectrum of transmon whose Hamiltonian is given as 


'^A = ^-Ej cos (C8) 

choosing to solve for the spectrum in the flux basis 
{|4>j)} where Qj = f we find 


As discussed in [27] He has a diagonal representation 
in terms of cosine basis. However, by choosing this basis 
H'^od remains as an A^-like term giving rise to intermode 
interaction. 

The solution to the above equation is a set of real 
eigenenergies and eigenmodes {hH„, 4'„((^j)|n € N°} 
Appendix C: CANONICAL QUANTIZATION OF where any operator in the transmon’s space has a spec- 
A CLOSED cQED SYSTEM representation over them as 


<P 

'2Cj 


— Ej cos ( 27r^ 


4>o /. 


(C9) 


We have all the tools to extend the classical variables 
into quantum operators by introducing a set of creation 
and annihilation operators as 

d{x, t) = X! (andn{x) + dld^(x)J (Cl) 


where 0{x,t) is any arbitrary classical field that we al¬ 
ready know its set of classical eigenmodes {On{x)} and 
eigenfrequencies {a;„} and 0{xH) represent its quantum 
analog. CpitOn) represents the appropriate normalization 
constant for each mode. The creation and annihilation 
operators obey the usual bosonic commutation relations 


dm d\P\ 

— iflSnm 

(C2) 

Qjfl , O-TTl] 

= 0 

(C3) 


= 0 

(C4) 


Remembering that {$„(a:)} represent Hermitian 
modes and thus real functions, we find the quantum op- 


Orit) =U{t)dT(0)U\t) 

=U{t) ( ^ (ml dT(0) |n) U\t) 

\m,n / ^ 

= ^ (m| Ot(0) \n) Pmn(t) 

m,n 


where {Pmn = |w) (n|} is a set of projection operators 
between states m and n , and 


{m\OT{0)\n) = / d(t)j'^rn{(l)j)OT 


4>j^ 


h d 
i d(t>j_ 


(Cll) 


We are now able to express and Qj in their spectral 
representation. Notice that since the potential is an even 
function of (pj, the eigenmodes are either even or odd 
functions of (pj, so the diagonal matrix elements are zero, 
since 
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{n\ |n) = / = 0 (C12) 

J '• - V- ^ 

Odd 

(n| Qj \n) = ^J d^j = Q (C13) 

'-^ 

Odd 

Therefore, we can express and Qj as 


= X! !”■) Pnin{t) 

m^n 

= J2 ("*i^./(o)i«) ( 

m<n 

Qjit) = X! \n) Pmn{t) 

m^n 

= ^ (m| Qj(0) \n) {Pmn{t) - Pn7n{t)j 

m<n 


(C14) 


(C15) 


where the second lines are written based on the observa¬ 
tion that by working in a flux basis, matrix elements of 
Qj and are purely real and imaginary respectively. 
Now that we know the spectrum of both the resonator 
and the qubit, we can easily write the interaction term 
as 


iQj 


dx 


p{xQ) 


6{x - xo) = 


c(x,xo) 

1 

^ ^ QJ,mniPmn Pnm) f ^ ^/(^o) 

m<.n,l ^ '' 


(C16) 


Defining the coupling intensity gmni as 


1 

2 

Wj ,mn )^l{xo) 


the interaction takes the form 


^ ^ ^9ninl{Pmn Pnrn){dl dj) 


d^rnnl — T 


/ hLOl \ 

\2cL ) 


m<n,l 


(C17) 


(CIS) 


Finally, up to a unitary transformation ai —>■ idi and 
Pmn —t iPmn for m < Ti, the Hamiltonian reads 


l~i — ^ ( dPljiPjiYi T ^ ( diOjid^Qn 


■Ha 




-b E ^9mnl (^^mn ^nm 


(C19) 


m<n,l 


Ui. 


Hamiltonian 

P = + X] diVn&ldn 



+ '^dgn{d- + d~'~){dn + al_) 


■Hi^t 


(C20) 


where we have used the shorthand notation a = Pqi, 
= Pio and wqi = fli — Hq. gmni is also reduced to 
g-a = 5oin given as 


%n = 7 



(jQj,oi)^n(a;o) 


(C21) 


In the main body of this paper we have only included 
the results for the case Xq = O.OIL. In Figs. |^and[^ we 
have considered the cases of xq = 0.50L and xq = 0.25L 
respectively. For xq = 0.50L we observe that all even 
numbered CC modes are unperturbed, while for odd 
numbered CC modes, both eigenmodes and eigenfrequen- 
cies are found to be less than cosine ones. The reason 
for invariance of even numbered modes is that originally 
qubit sits on a local minimum of the photonic energy den¬ 
sity and therefore does not interact with these modes. In 
addition, regardless of the value for Xs-, ah odd numbered 
coupling strengths are zero due to high symmetry of this 
point. For Xs = Oj the even modes follow an envelope 
that goes like and as Xs grows they are suppressed 
and fall below this envelope. 

This behavior is not specific to xq = 0.50L. Generally, 
if the qubit is placed at Xq = G N then there is 

a periodicity in the mode structure such that every n 
modes remain unperturbed. This is for example observed 
in Fig. where a;o = 0.25L and thus modes indexed as 
4n — 2, n S N are unchanged. In this case, depending on 
the value of Xs the CC coupling strengths can be either 
below or above the solutions for Xs = 0. 


Appendix D: TRK SUM RULES FOR A 
TRANSMON QUBIT 


In this Appendix, first we are after finding a general 
sum rule in quantum mechanics and then we apply the 
results to calculate upper bounds for matrix elements of 
charge and flux operator i.e. (m| Qj \n) and (m| <!>,/ |n) 
for the case of a transmon qubit. Assume a Hamiltonian 
% where its eigenmodes and eigenenergies are known as 
{|n) ,£'„|n € N*^}. Consider an arbitrary Hermitian op¬ 
erator A = where we define successive commutation 
of 'H and A as 


Notice that by truncating transmon’s space into its 
first two levels, we are able to recover a multimode Rabi 


d(^) _ 
- 




(Dl) 
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FIG. 6. Closed-boundary CC modes for xo = 0.51/ (a)- First 
20 Eigenfrequencies for Xs =0.1 b-e) Normalized energy den¬ 
sity of the first 4 modes for Xs = 0.1. The black curve shows 
cosine modes while the red ones represent CC modes. The 
blue star shows where the qubit is connected, f) First 20 
coupling strengths for various values of Xa- 


FIG. 7. Closed-boundary CC modes for xc = 0.25L (a)- First 
20 Eigenfrequencies b-e) Normalized energy density of the 
first 4 modes for Xs = 0.1. The black curve shows cosine 
modes while the red ones represent CC modes. The blue 
star shows where the qubit is connected, f) First 20 coupling 
strengths for various values of Xs- 


By this definition we find that for any two arbitrary 
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eigenstates \m) and \n) 

{m\ |n) ={Em - En) (m| \n) 

= {Em - En)’" {m\A\n) 
Using the above identity we find that 


Noticing that = 2e brings us the sum rule for Qj as 


(fc) 


A 


\m) = {m\A \m) 

EU><"I 

— (m|C^^ A\m) 

El«>("l 


(D2) 


(D3) 


= Y,iEn-Em)[l-{-l)'^] {m\A\n) 


Now, consider the Hamiltonian for a transmon qubit 


00 


Y^^iEn- Eo) 


n>0 


^ / Q-TT 

\Qj\n) = (2e)^£’j (0|cos f —$j ) |0) 

< {2efEj 


(DIO) 

Again, due to positivity of terms on the L.H.S we hnd 
2e'^ET 2Ej 


IQj.oil < 


Ui - Uo ^ VSUjUc - Ec 


(Dll) 


Appendix E: GENERALIZATION TO AN OPEN 
cQED SYSTEM 

1. Lagrangian and Modified Eigenmodes 




(D4) 


where C/(4>j) = — Ej cos . Applying the result 

found in ID3I we can write 


(D5) 


[$J, [H, $j]] |0) = ^ 2 (U„ - Uo) (0| l>j \n) 


/^(l) 


n>0 


where |0) represents the ground state. The L.H.S can be 
calculated explicitly as ^ which leads to the sum rule 
for as 

5]2(if„-ifo) (0|4>j|n) ' = ^ (D6) 

Noticing that all terms on the L.H.S are positive, we can 
find an upper bound for 4)j oi as 


I^,7,011^ < 


^2 


Ec 


n 

2Cj{Ei—Eq) A^EjEq — Ec V® 


(D7) 


where we have defined the charging energy Ec = 

In a similar manner, it is possible to obtain a sum rule 
for Qj as 


[Qj, [n, Qj]] ]0) = Y^ - ^o) Qj ]n) 


^(1) 


n>0 


(D8) 


The L.H.S can be explicitly found as 


[Qj, [E, Qj]] = E 






V 00 / 


, 27r » 
Ej cos — 


(D9) 


The results from the previous section make it very easy 
to find the Lagrangian and hence the dynamics for the 
open case where now the end capacitors Cr and Cl have 
finite values as shown in Fig |lb| Here, we have a fi¬ 
nite length resonator which is capacitively coupled to two 
other microwave resonators at each end. Assuming that 
the transmon qubit is connected to the resonator at some 
arbitrary point x = Xq, the Lagrangian for this system 
can be written as 


C = 

+ 

+ 

+ 

+ 

+ 

+ 




[L 

/ dx 

lo+ 

pLA-Ln 

Il+ 

r 

J-L-Ll 


I 

2^dt’ 2Edx’ 


dx 


dx 


0^77x2 1.0^77x2 


2'< dt 

'1 

Asf 


,x2_ 

2E dx ’ 

^ 2r dx ^ 
2 


Ui($i.(0-,t)-$(0+,t)) 
Ufl($fl(L+,t)-$(L-,t)) 
Cg (^j(t) - ^(a:o,i)) 


(El) 


We have already learned how the coupling intensity 
depends on the Hermition eigenmodes and eigenfrequen- 
cies of the resonator as well as the dipole moment of the 
transmon. Here we have the same situation except that 
due to the opening introduced by the finite end capaci¬ 
tors Cr and Cr, we need to find the modified Hermitian 
modes of the open system. Therefore, let’s for the mo¬ 
ment forget about Lagrangian of the transmon and its 
coupling to the resonator and focus on the modification 
introduced by one of the end capacitors, for instance Cr- 
The trick is that we can write this contribution as sum 
of three separate terms 
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\cl (‘i>L(0-,i) - 4>(0+,t))' = 

+ iCi$(0+,t)2 

-CL$(0+,i)$L(0-,t) 

(E2) 

Notice that only the term ,t) is re¬ 

sponsible for coupling of the resonator to the bath and 
the other two can be considered as a modification on top 
of the closed case. Applying the same method for the 
right capacitor, we can define the modified Lagrangian 
for the left and right baths 






dx 


lL + 


1 


2 '’* at 


’ 2r dx ’ 


-CR^RiL+,t)^ 


dx 


1 

r^si 


1 ,d^T 


^ ^ ^ 2V dx ’ 


(E3) 


(E4) 




In addition, the interaction Lagrangian is found as 

^C,LR = —C'L‘h( 0 ^, <)<i)L( 0 “, t) — Cr^{L~ ,t)i^R{L'^ ,t) 

(E5) 

which is also equal to the interaction Hamiltonian since 
by going from Lagrangian to Hamiltonian capacitive con¬ 
tributions (kinetic contributions) don’t change sign. The 
idea is to find the Hermition modes governed only by 
each of these uncoupled modified contributions and fi¬ 
nally write the interaction in terms of Hermitian modes 
of each subsystem. 

Up to here, we haven’t considered the effect of trans- 
mon on capacitance per length as we found in |B19| It 
is not necessary to go over the derivation again, since 
these two effects, modification due to opening and due 
to transmon, are completely independent. By consider¬ 
ing the inhomogeneity introduced by the transmon we 
have 



1 / n/54>,o 1 ,94),o 

/ dx 

/o+ 





where Cop(x,xo) is given as 

Cop(x, xo) = c + Cs5(x - Xo) 


+ CrS(x-L-) + ClS(x-0+) 


(E 6 ) 


(E7) 


The new delta functions in Cop(x) are only important 
at the boundaries, which can be found by integrating the 
equation along an infinitesimal interval that includes the 
delta functions. In order to find the modes, we need to 
solve 


d^^n(x) 

dx^ 


+ lcLJ^^n(x) =0, X ^ Xo 


with boundary conditions given as 
= lCRu;^i>n(L-) 




x—L 


^n(a:) 


d_ 
dx 

d^n{x) 


= -lCLUjl^rt{0+) 


(E 8 ) 

(E9) 

(EIO) 


x= 0 + 


dx 


d^n{x) 


dx 


^n{xQ ) — $n(a^Q ) 


ICstul^nixo) = 0 (Ell) 


(E12) 


Defining unitless parameters xr,l = - 3 ^ as we did 
for Xst we find eigenfrequencies satisfy a transcendental 
equation as 


+ (1 - XRXhiKLf) sin (fc„L) 

+ {XR + Xl) knL cos (fc„L) 

-I- Xsk„L cos (knXo) cos (fc„(L - Xq)) 

- XRXsiknL)"^ cos (knXo) sin {k„{L - xq)) 

- XhXsiknLY sin {k^xo) cos (/c„(L - a;o)) 

+ XRXLXs{knLf‘ sin (fc„xo) sin {kn{L - xo)) = 0 


(E13) 


and the real-space representation of these modes read 

(E14) 


$„(x) oc 


^n{x) 0 < X < Xo 

^n(x) Xo < X < L 


where 4’<(a;) and 4>>(a;) are found as 
l><(a:) = [cos(fc„(L-a;o)) - XRknL sin {kn{L - xo))] 

(E15) 


X [cos (knx) - XLknL sin (knx)] 


$>(a;) = [cos(fc„xo) - Xi/c^Lsin (fc„a;o)] 

X [cos {kn{L -x)) - XRknL sin {kn{L - x))] 

(E16) 

The normalization constant will be set by the orthogo¬ 
nality conditions 


^^ Cop(x, Xo) 


dx^^^^^^ i9^n(a;) 


(E17) 


dx 


dx 


1 
2 

^nLSmn 


- O (*- + L XR^nt{L-)^n{L-) + (0+)l>„ (0+) 

(E18) 
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2. Canonical Quantization 

Following the same quantization procedure as for the 
closed case we can write the field operators in terms of 
the eigenmodes and eigenfrequencies of each part of the 
circuit as 


71 


n 


2uJr, cL 


[an{t) + a\it)) l>„ (x) (E19) 


(E20) 

(E21) 

where we have also considered some finite length for 
the left and right resonators as well to keep the normal¬ 
ization constants meaningful. Now, we can easily de¬ 
rive an expression for resonator-bath coupling in terms 
of modes of each part. Consider coupling to the left bath 
for the moment 


The expression for gn is the same as in |C17| but with 
the new set of Hermitian modes satisfying the open¬ 
boundary conditions discussed before. The interaction 
Hamiltonian then is found as 

^ ^ b^9mnl{Pmn Pnrn){bl dj) (E28) 

m<n,Z 


Gathering all different contributions together and mov¬ 
ing to a new frame where Pmn —t iPmn for m < n, 
an —>■ ittn and hn^L/R ibn^L/Rj the Hamiltonian in 
its 2”^^ quantized form reads 


H = E hj^riPrin H” E hWndidn -I- Y f^^n,sb\^s^n,S 

n n n,S={L,R} 


nA 


Uc 




^ ^ ^Qmnl ^ 


Pmn Pnm 1 ( 


m<.n,l 


y-i 


“ 1 “ ^ ^ ^Pmn,S (^m “t“ ^In) (j^n,S P ^n,S^ 


m,n,S—{L,R} 


y-CE 


(E29) 


7icL = -GLd(0+,t)dL(0-,t) (E22) 

ncR = -CR^L+,t)^R{L-,t) (E23) 

Considering that the time-dynamics of annihiliation 
and creation operators up to here are only governed by 
the free Lagrangian of each part, we have 

t) =Y (^ 2uj^ cl ) (t)) ^n(a;) 

= “*E(^) (dnit) - di{t)) ^n{x) 

(E24) 


And we have the same type of expression for the end res¬ 
onators as well. The interaction Hamiltonian then reads 

Pc,LR — 'y ^ df^nin^R (dm ^m) (bn,R ^n,Rj 


m,n 

^ ^ ^(^mn,L (^m ^m) ^n,L^ 


(E25) 


m^n 

where we find Pmn,R and ^mn,L as 


l3mn,R = -^^^=UjkujlR^m{L-)^nAL'^) (E26) 

An,L = (E27) 


Appendix F: LARGE TRANSMON LIMIT 


It is very insightful to see how the previous results 
change if we assume transmon’s length d is comparable 
with resonator’s wavelength L. In such a case we assume 
that the coupling is not local and spreads over whole 
length of transmon with mutual capacitance per length 
Cg. Eollowing the same discrete to continuous approach 
one finds the Lagrangian as 


C=-{Cj + Cgd)A-Uj{<^j) 


dx 


1. . 

-(c + c,M».-o))(^) -y(^) 


dx CgTTd{x,Xo)^J-^ 


(FI) 


where TTd{x, xq) is a unit pulse of width d that is defined 
in terms of Heaviside function 0{x) as Trd{x,xo) = 9{x — 
xq + d/2) — 9{x — Xq — d/2). The Euler-Lagrange E.O.M 
then read 


{Cj + Cgd)^J + 


dUjj^j) 

d^.j 



5^$ 

dxCgTid{.x,Xo)-^ 


(F2) 

-l{c + CgTTdix, Xo)) = -IcgTTdix, a;o)$j (F3) 

Comparing this to the Euler-Lagrange E.O.M drived 
earlier in appendix B.2 It is clear that the structure of 
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the equations has remained the same while Cg5{x — Xg) 
has been replaced with Cg7r(x, xg) ■ The Hamiltonian can 
be found through the usual Legendre transformation as 


n = 


Q 


2 

J 


2C 


mod 

J 




/.L 

/ dx 

Jo 

p^{x,t) ^ 1 


2cd{x,Xg) ' 21 

\ i9a; / 

Qj 

dxCgTTd{x,Xg) 

0 

p{x,t) 

CJiod J 

Cdix,Xo) 



(F4) 


where Cd{x,Xo) is the modified capacitance per length 
and reads 


Cd{x,Xo) 


c + CgTrd{x,Xo)® 


d 


(F5) 


1. Classical Lagrangian 

Following the usual canonical quantization scheme, 
first we have to find the classical Lagranian. We already 
know the equations of motion for the EM fields to be the 
Maxwell’s equations as 

V.E(r,t) = (Gl) 

eo 

V.B(r,t) = 0 (G2) 

VxE(r.,) = -^ (G3) 

V X B(r,0 = (G4) 

where p(r,t) and J(r,t) are scalar charge density and 

vector current density and are given as 

P(r,t) = J^g„d(^Ur-r„(t)) (G5) 

n 

J(r,t) = (r-r„(t)) (G6) 


The ©-notation represents series combination of two ca¬ 
pacitors. Surprisingly, we observe that when the dimen¬ 
sion of transmon is taken into account, the modification 
is mutual and transmon’s spectrum is also influenced by 
the coupling such that reads 


^mod ^^ 


cCg7rrf(a:,Xo) 
/o c-|-Cg7rd(a:,a;o) 


Cg-Rj,(x,Xo)®C 


= Cj + {cg®c)d 


(F6) 


The remaining equation of motion is a Newton equa¬ 
tion regarding the mechanical motion of the electron 
which reads 


TO„r„(t) = Qn [E(r„(t),t) -k r„(t) x B(r„(t),t)] (G7) 

'-V-' 

Lorentz Force 

Based on |G2| and |G3| we are able to express the phys¬ 
ical fields E(r,t) and B(r, t) in terms of scalar potential 
V{r,t) and vector potential A(r,t) up to a gauge degree 
of freedom as 


The results above are general such that one can re¬ 
place the rectangular pulse CgTrd{x, ccO) in capacitance per 
length by any smooth capacitance per length Cg(x,Xo) 
and the form of ?{, Cd(x, Xg) and remain the same. 


Appendix G: HAMILTONIAN AND MODIFIED 
EIGENMODES OE A CLOSED CAVITY-QED 
SYSTEM 


In this Appendix, first we derive the classical Hamil¬ 
tonian for a general system containing finite number of 
point charges interacting with the EM field inside a closed 
cavity. This is achieved by expressing the Maxwell’s and 
Newton’s equations of motion in a Lagrangian formalism 
and then a Legendre transformation to find the Hamil¬ 
tonian. This model is then reduced to describe a one¬ 
dimensional cavity shown in Fig. In order to empha¬ 
size on the resemblance to the cQED results we found 
earlier, it is assumed that there is only a single electron 
at Re(t), while all other electronic or nuclear degrees of 
freedom are frozen at Rq. 


E(r,t) = -^^^-VE(r,t)) (G8) 

B(r,t) = V X A(r,t) (G9) 


It is possible to write a Lagrangian that produces all 
previous equations of motion G1||G4 and CT as a result of 
the variational principle SjC = 0. This Lagrangian reads 




m„r^ 


d^r ^£o(dtA + VVf - ^(V X Af (GIO) 

_2 ZiiQ 


+ J d^r[J.A-pV] 


In order to proceed further, we need to fix the gauge. 
Choosing to work in Coulomb gauge defined as V.A = 0 
and using |G1| we find that the scalar potential V (r, t) 
satisfies a Poisson equation as 


VV(r,t) =(Gll) 
eo 
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Having the charge density as G5 we can easily solve 
this equation to obtain 


substituting |G16 and G17 into the expression for 
Hamiltonian we hnd 




E 


q-n 

47reo|r - r„(t)| 


(G12) 


Furthermore, this choice of gauge helps to simplify the 
Lagrangian since due to the divergence theorem 

J d^rdtA.WV = J d^r\7.{VdtA) - j d^rVdt (V.A) 

0 

= ^ dS. {VdtA) 

(G13) 


which means this term only contributes at the bound¬ 
aries. Boundary terms do not affect equations of motion 
inside the cavity, however their existence are necessary to 
ensure the correct boundary conditions, i.e. continuity of 
parallel electric held and perpendicular magnetic held at 
the interface of cavity with the outside environment. As 
far as we hx these conditions properly, we can remove all 
surface terms form the Lagrangian. In a similar manner 


= E + E 


-eo(i9tA(r,t))2 + (V x A(r,t)) 


2/xo 


(G19) 


By replacing 9tA(r,t) and r„(t) in terms of conjugate 
momenta r[(r, t) and Pnit) respectively, the Hamiltonian 
can be rewritten as 




+ / d\ 


n2(r,t) , (VxA(r,t))^ 


2en 


2^0 


(G20) 


J dVieo(VH)2 = ieo f dS. (VW) - J d^r^coVV^V 

= ieo^ds.(yvy) + J d\^pv 

(G14) 

Finally, by putting everything together and neglecting 
surface terms, we hnd the simplihed Lagrangian as 




+ J d\ ieo {dtAf - ^ (V X A)^ 

-h f d^rJ.A 


(G15) 


2. Classical Hamiltonian 


The hrst step is hnd the conjugate momenta as 


which can be written in a more instructive way as 


n 

n 



Ha 



n2(r,t) 

1 (VxA(r,t)f 

2eo 

2^0 





^ (G21) 

+ j d^r^^^A^{r,t)5^^\r-Yn) 

'■- - -V-' 

_E i!Lp^_A(r„,t) 

'-V-" 


d/L 

Pn = =mnrn + qnA{rn{t),t) (G16) 

dC 

^ = eo9t A(r, t) (G17) 

i9A 

Then, the Hamiltonian is calculated via a Legendre 
transformation of the Lagrangian as 


'H = ^p„.r„-t- d^rn{r,t).dtA{r,t) - £ (G18) 


This is the most general form of the classical Hamil¬ 
tonian of a hnite number of charges interacting with the 
EM held inside a closed cavity. In what follows, we make 
a few assumptions to reduce this model for the system 
shown in Fig. First of all, we assume that the wave¬ 
length of EM held is much larger than atomic scale 
such that we can apply 0*^-order dipole approximation 
A(Re(t),t) « A(Ro,t) in both and T-Lint- Finally, 
by considering that there is only a single charge located 
at Re(t) and the nucleus is hxed at Ro the Hamiltonian 
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is reduced to 




■Ha 


+ / d^r 


nHv,t) , (VxA(r,t))' 


2en 


2^0 


He 


+ [ d\^A\r,t)6^^\r-Ko) 

J ^rUe 




-Pe.A(Ro,t) 

TOe 


Hi. 


(G22) 


Assuming a closed cavity case and remembering that 
E|| and Bj^ are continuous across the cavity, the bound¬ 
ary conditions read 


ny X za;A(r, w)-b 
nj_. X A(r,w)^ 


= 0 


= 0 
B 


(G28) 

(G29) 


where nj_ and ny represent perpendicular and parallel 
unit vectors on the boundaries of the cavity and V is the 
time-Fourier transform of the scalar potential. Equation 
|G27| with the boundary conditions above provide a dis¬ 
crete set of modes due to finite volume of the cavity. For 
notation simplicity, we label the eigenfrequencies as uj\ 
and the modes as AA(r) = A(r,a;A) , while in reality A 
denotes multiple sets of discrete numbers each for a sep¬ 
arate dimension of the cavity. These modes satisfy the 
general orthogonality relation 


3. Modified Cavity Eigenmodes and 
Eigenfrequencies 

Having found the modification introduced by the A^ 
term in the previous section, we can now calculatae the 
effect it has on the structure of the modes. Specifically, 
we are after eigenmodes of the modified Hamiltonian for 
the cavity given as = He + H^°'^. This can be 

done by finding the Hamiltonian E.O.M for the conjugate 
fields as 


^A(r,t) = —n(r,t) 

Ot Co 

4n(r,t) = —X (V X A(r,t)) 

Ot fJ,Q 

- —A{r,t)S^^\r - Ro) 
trie 

By combining these two equations and applying the 
gauge condition V.A = 0 we find 

- Moeo A(r, t) = A(r, (r - Rq) 

(G25) 

which is a wave equation with an extra term on the R.H.S 
due to A^ modification. The coefficient can be ex- 

rrie 

pressed in terms of fine structure constant a and Bohr’s 
radius Oq as dTra^og. By doing a Eourier transform 

1 

A(r,t) = — dujA{r,uj)e (G26) 

27^ 7-00 

we can easily separate the time and spacial dependences 
to obtain 

^ — 47ra^ao(J^^^(r — Rq)^ A(r, w) = 0 

(G27) 


(G23) 

(G24) 



V<5aa' 


(G30) 


where we have set the normalization such that the modes 
are dimensionless. Another orthogonality relation can be 
found in terms of VAA(r) as 


+ / d^rVAA(r).VAA'(r) 


dS. 


AA(r).VAA'(i’) + AA'(r).VAA(r) 


(G31) 


47ra^aoAA(Ro).AA'(Ro) = kxkyV6xx> 


Up to this point, we have considered the mode struc¬ 
ture for a general cavity with any arbitrary geometry. 
In order to demonstrate the connection to the results 
for a one dimensional cQED system, we have to make 
a few assumptions about the geometry of the cavity. 
We assume that the cavity’s length is much larger than 
the diameter of its cross section, i.e. L ^ ^/S. By 
considering variation of the eigenmodes only along this 
dimension we can write A(r, t) = \izA(xH) and thus 
B(r,t) = —\iydxA{xH)- The Hamiltonian is then re¬ 
duced to 




Ha 


[ d?s [ dx 

lP{x,t) 1 {dxA{x,t)f 

J Jo 

2eo 2fiQ 


He 


+ [ d^s [ dx— —A^(a;o, — So)i5(a; — a;o) 

J Jo 2TOe 




- 

trie 


Hi. 


(G32) 
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Following the same procedure, we can find a modified 
wave equation as a result of 

+ (“)' - - -.)) = 0 (G33) 

Assuming that the atom is fixed at point xo,W only 
affects the bouundary condtion for the zero frequency 
mode which we are not interested in. Therefore, by ap¬ 
plying the boundary conditions 


A(x, oj) 


x—0,L 


= 0 


(G34) 


we find the normalized eigenfrequencies to satisfy a tran¬ 
scendental equation as 


. sm(knXo)sm(kn(L - xo)) „ 

sm (knL)+Xc - r~F -=0 (G35) 

k„L 


where we have defined the unitless parameter Xc = 
47ra^oo§. The real-space representation of the eigen- 
modes read 


An{x) OC 


sin (fc„(L — xq)) sin (fc„a;) 0 < a; < Xg 

sin (knXo) sin (fc„(T — x)) xq < x < L 

(G36) 


Eventually, one can show that these eigenmodes satisfy 
the orthogonality relations: 


dxAm{x)An{x) = L5„ 


(G37) 


/o 


r)A r)A V ~ 

+ ^Am(a;o)A„(xo) = kmknL5„ 


/o 


dx dx 


(G38) 


Note that only the ratio ^ is determined by the ge¬ 
ometry of the cavity, while the pre-factor ~ 

3.54 X is a universal length scale. This implies 

that the modification is only visible when £ is around 
the same order as Ana^ao. 


4. Canonical Quantization 


Having found the proper set of eigenmodes and eigen¬ 
frequencies that diagonalizes the classical Hamiltonian 
for a one-dimensional closed cavity-QED system, we can 
move forward and extend the classical variables into 
quantum operators by introducing the necessary commu¬ 
tation relation between conjugate pairs. Let’s consider 
the conjugate fields for the cavity first. We can expand 


these fields in terms of the proper eigenmodes as 

I 

= + (G39) 

n(x,t) =(a„ - oj;) i„(x) (G40) 

where a„ and ajj are annihilation and creation oper¬ 
ators for each mode. By inserting the above equations 
and using the orthogonality conditions |G37| and |G38[ 
'fl'fnod _ _|_ ymod become diagonal as 

rj^rnod ^ ^ anal) = ^ tu^nolan + COUSt. 

n n 

(G41) 


The next step is to obtain the spectrum of Ha by solv¬ 
ing a Schrodinger equation in real-space basis as 

- eE(re)^ 4'n(re) = (G42) 

where we have denoted the eigenmodes and eigenenergies 
by {4'„(re),E„ = hi}n\n G N°}. Ha can be decomposed 
as 


Ha = ^ hVtnPn 


(G43) 


Pe also has a spectral decomposition over this basis. 
Since the Goulomb potential V{re) is an even function 
of Tg, as we explained in the case of charge qubit only 
diagonal matrix elements of pg are nonzero and we can 
write 

Pe = ^ ^ (w| Pg |n) Pmn (G44) 

m^n 

where matrix elements Pe,mn can be calculated as 

/ t 

d^rg4'm(re)^V4'„(rg) (G45) 

By working in a basis where 'I'„(rg) are real functions, 
the dipole matrix elements are purely imaginary which 
allows us to write 

Pe = 'y ^ (w| Pe 1*^) {Pmn Pnm) (G46) 

m>n 

Finally, for sake of resemblance to the cQED results we 
move to a new frame Pmn iPmn for m < n to obtain 
the Hamiltonian as 

H — y ^ kPlyiPrfiyi -f y ^ flCOndmdn 


-Ha 




T y ^ ^9mnl (^Pmn H Pnm^ 4” ^ 


Tn>n,l 


T-iin 


where the coupling strength gmni reads 


k>9mnl — (^Pe,mn-Cl2) . g, ^ 

TOg \2eoU!lbL 


(G47) 


Mxo) (G48) 
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